Optimal fidelity of teleportation of coherent states and entanglement 
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We study the Braunstein-Kimble protocol for the continuous variable teleportation of a coherent 
state. We determine lower and upper bounds for the optimal fidelity of teleportation, maximized 
over all local Gaussian operations for a given entanglement of the two-mode Gaussian state shared 
by the sender (Alice) and the receiver (Bob). We also determine the optimal local transformations 
at Alice and Bob sites and the corresponding maximum fidelity when one restricts to local trace- 
preserving Gaussian completely positive maps. 
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I. INTRODUCTION 

Quantum teleportation [l( is the transfer of an un- 
known quantum state from a sender (Alice) to a receiver 
(Bob) by means of the entanglement shared by the two 
parties and appropriate classical communication. Bob 
recovers an exact copy of the state teleported to him by 
Alice only if the quantum channel is the ideal maximally 
entangled state, which however, in the case of continu- 
ous variables (CV), is an unphysical infinitely squeezed 
state Nevertheless, by considering the finite quan- 
tum correlations between the quadratures in a two-mode 
squeezed state, Braunstein and Kimble [1J proposed a re- 
alistic protocol employing a beam splitter and homodyne 
measurements, which approaches perfect teleportation in 
the limit of infinite degree of squeezing. This protocol 
(and its various extensions) has been then implemented 
by various groups [H . The success of a teleportation ex- 
periment is quantified by the fidelity which, in the case of 
a pure input state \ipi n ), is given by T = (ipin\pout\^in) , 
where p ou t denotes the output state of the protocol, and 
coincides with the probability of finding the input state 
\ipi n ) at the output. 

The relation between the fidelity of CV teleportation 
T and the entanglement shared by the distant parties 
is nontrivial and it has been investigated in various pa- 
pers H, @. fn fact, entanglement is the key resource 
that allows to beat any classical strategy for transmit- 
ting quantum states. One can say that genuine quan- 
tum teleportation has been performed only if T > J- C i, 
where F c i is the classical fidelity threshold achievable by 
two cheating parties who can perform arbitrary local op- 
erations and classical communication (LOCC) but are 
not able to share entanglement, nor to directly trans- 
mit quantum systems [5]. This means that T > T c i 
is a sufficient condition for entanglement between Alice 
and Bob; it is not a necessary condition because one can 
have lower-than-classical fidelities even using an entan- 
gled state. This is due to the fact that the Braunstein- 
Kimble protocol chooses specific combinations of quadra- 
tures, and therefore T , differently from entanglement, is 
not a local symplectic invariant: if the protocol employs 
combinations of quadratures which are inappropriate for 
the given shared entangled state, the resulting T is very 



low. The difference between of teleportation fidelity and 
entanglement opens two problems: i) the determination 
of the classical fidelity threshold T c i for a given class of 
CV input states; ii) the optimization of the teleportation 
fidelity, for the chosen class of input states, over all pos- 
sible LOCC strategies. Determining T c i is a nontrivial 
quantum estimation problem which has been solved only 
for few classes of states: i) input coherent states with 
completely unknown amplitude {T c i = 1/2) [f|; ii) pure 
squeezed states with zero displacement and completely 
unknown degree of squeezing at a given phase (see Ref. p} 
which obtained T c i = 0.81 by considering LOCC strate- 
gies in which one is only allowed to prepare squeezed ther- 
mal states); iii) pure squeezed states with completely un- 
known displacement and orientation in phase space but 
fixed degree of squeezing s {T c \ — y/s/(l+s)) @|. Here we 
restrict to input coherent states, which represent the ba- 
sic resource for many quantum communication schemes 

The improvement of the teleportation of coherent 
states by means of local operations and its relation with 
the entanglement of the shared entangled state has been 
already discussed in a number of papers [f| [II], [H, [Uj] . 
Ref. [Il( showed that in some cases the fidelity of tele- 
portation may be im prov ed by local squeezing transfor- 
mations, while Ref. [l2] showed that in the case of a 
shared asymmetric mixed entangled resource, teleporta- 
tion fidelity can be improved even by a local noisy oper- 
ation. Ref. [l3( then considered the class of local trace- 
preserving Gaussian completely positive (TGCP) maps 
(those performed by first adding ancillary systems in 
Gaussian states, then performing unitary Gaussian trans- 
formations on the whole system, and finally discarding 
the ancillas), and maximized the fidelity over this class 
of operations. 

Ref. [l3|] confirmed that the best local TGCP map 
maybe a noisy one, i.e., that teleportation fidelity can 
be increased even by decreasing the entanglement and in- 
creasing the noise of the shared entangled state. Ref. 
however, did not discuss the relationship between entan- 
glement and the optimal fidelity T op t- Ref. § instead 
found this relationship, but only for a subclass of sym- 
metric Gaussian entangled state shared by Alice and Bob: 
for this class it is T op t = (1 + v) ~ , where v is the low- 
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est symplectic eigenvalue of the partial transposed (PT) 
state. The parameter v provides a quantitative charac- 
terization of CV entanglement, because the logarithmic 
negativity E^f is related to v by E^f — max[0, — Inv] [ijj]. 

The problem of finding the optimal LOCC strategy is 
non-trivial because the set of operations that Alice and 
Bob can adopt is very large. In fact, apart from local 
TGCP maps, they can adopt two further options: i) use 
non-trace preserving Gaussian operations in which some 
ancillary mode is subject to Gaussian measurement, i.e., 
projected onto a Gaussian state, rather than discarded 
[24|; ii) use local non-Gaussian operations (either with 
measurement on ancillas or not), i.e., those involving in- 
teractions which are non-quadratic in the canonical co- 
ordinates. The first class of maps, together with TGCP 
maps, forms the most general class of Gaussian com- 
pletely positive (GCP) operations, capable of preserv- 
ing the Gaussian nature of the state shared by Alice and 
Bob. Non-Gaussian operations instead will transform the 
initial Gaussian bipartite state of Alice and Bob into a 
non-Gaussian one, and they can also increase the fidelity 
of teleportation in some cases [2|| . 

In this paper we generalize in various directions the 
results of Refs. [H, [lj|. We show that if Alice and Bob 
share a bipartite Gaussian state with a given v and one 
restricts to local GCP maps which preserve such a Gaus- 
sian nature, the optimized fidelity always satisfies 



where 



1 



l + 3f 
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< 



1 



1 



(1) 



We also show that the upper bound is reached iff Alice 
and Bob share a symmetric entangled state. Moreover 
we determine the optimal local transformations at Alice 
and Bob sites and the corresponding value of J- op t as a 
function of the symplectic invariants of the shared CV 
entangled state when one restricts to local TGCP maps. 

The paper is organized as follows. In Sec. II we provide 
the basic definitions of the problem, in Sec. Ill we prove 
and discuss the lower and upper bounds for the optimal 
teleportation fidelity for a given shared entanglement. In 
Sec. IV we discuss the properties of the optimal local map 
and derive its explicit form in the case of TGCP maps. 
Sec. V is for concluding remarks. 



II. DEFINITION OF THE PROBLEM 

The protocol for a perfect CV quantum teleportation 
based on ideal Einstein-Podolski-Rosen (EPR) correla- 
tions has been introduced by Vaidman 0, and then 
adapted to the finite correlations of a two-mode squeezed 
states by Braunstein and Kimble [|| . The idea can be eas- 
ily shown in the Heisenberg picture. We consider two CV 
systems, each described by a pair of conjugate dimension- 
less quadratures £k and pt (k — a,b). Introducing the 
vector £ T = (x a ,Pa,Xb,Pb), we can write the canonical 
commutation relations as [£/,£ m ] = iJim (l,m = 1, ...4), 



J=J®J, J = 



1 

-1 



(2) 



and ® denotes the usual direct sum operation. 

Alice and Bob share a bipartite CV state with EPR- 
like correlations, i.e., a state which can be considered 
as an approximate simultaneous eigenstate of the com- 
binations of quadratures x a + Xb and p a — pb, so that 
the variances of these two combinations are both close to 
zero. Alice also possesses an unknown input state with 
quadratures Xi n ,Pi n which she wants to teleport to Bob. 
Alice mixes the input mode with her part of the entangled 
state via a balanced beam splitter and she carries out a 
homodyne detection on each output mode, thereby mea- 
suring two commuting quadratures x+ = (x a + 5;, n )/v2 
and p- = {pin —p a )/\2. After receiving the measured 
values x + and p_ from Alice, Bob uses this transmitted 
classical information to perform a suitable conditional 
displacement on his own mode, Xb — ► x' b = Xb 

Pb — >p' b = Pb + V%p- ■ 

If we assume ideal homodyne detectors on Alice site, 
and that the shared bipartite state is undisplaced (i.e., all 
mean values of Alice and Bob quadratures vanish), the 
EPR-like correlations x a ~ —Xf, and p a ~ pb, together 



with Bob displacements, imply x' b 



and p' b ~ pi. 



i.e., Bob mode is described by a pair of conjugate vari- 
ables very close to those of the input mode. In the 
Schrodinger picture, this is equivalent to teleport the in- 
put state to Bob with a fidelity very close to one. 

We restrict to the case when the state shared by 
Alice and Bob p a b is Gaussian, where a compact ex- 
pression of the resulting fidelity of teleportation can 
be derived [H, EE B3- A bipartite CV state p a b is 
Gaussian if its Wigner characteristic function <I?afc(£) = 



Tr[p ai ,exp(-i£ i ■ £)] (where £ 7 = (x a ,p a , x bl Pb) is the 
vector of phase-space variables corresponding to £ T ), is 
Gaussian, i.e., <& a b(l) = exp(—£ r V£/4: + id?£). We have 
assumed that Alice and Bob share a zero-displacement 
state, implying d — 0. Therefore p a b is fully character- 
ized only by its correlation matrix (CM) V, whose generic 
element is defined as VJ m = (A£;A£ m + A£ m A£j) where 
A£/ = — (£/). The CM satisfies the uncertainty princi- 
ple V + ij > [13, and can always be put in the block 
form 



(3) 



where A, B, and C are 2x2 real matrices. Using charac- 
teristic functions, it is straightforward to prove [H, [l5[ 
that, if the input state is a single-mode Gaussian state 
with CM Vi n , the fidelity of teleportation is given by 

(4) 



^/det (2V m + N) ■ 



where 



N = ZAZ + ZC + C Z + B, 



(5) 
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with Z = diag(l, -1) [l8|. The 2x2 matrix N is semi- 
positive definite, N > 0, it describes the noise added to 
the teleported state, and it is equal to zero only in the 
ideal situation of perfect EPR correlations between Al- 
ice and Bob. As discussed in the introduction, we shall 
restrict to the case of input coherent states, Vi n — I, so 
that Eq. (TJ| reduces to 



T = 



V4 + 2TrA + det N 



(6) 



The problem afforded in this paper, i.e., the maximiza- 
tion of the teleportation fidelity over all possible Gaussian 
LOCC strategies for a given Alice-Bob entanglement, 
therefore means to determine the optimal local transfor- 
mation of matrices A, B and C which makes N as small 
as possible. 

As showed in [16j |, using an unbalanced beam splitter 
is equivalent, for the teleportation protocol, to a squeez- 
ing operation by Alice. Therefore the optimization over 
all Alice and Bob local operations includes also any even- 
tual modification of the beam splitter used for the joint 
homodyne measurement. 



III. UPPER AND LOWER BOUNDS FOR THE 
FIDELITY OF TELEPORTATION 

In this section we prove Eq. (fTJ), i.e., the upper and 
lower bounds for the fidelity of teleportation for input co- 
herent states. An important preliminary result enabling 
us to derive the two bounds is the fact that the optimal 
noise matrix N is very simple: in fact, the maximum 
teleportation fidelity is obtained when N is proportional 
to the 2x2 identity matrix /. More precisely, we have 
the following 

Lemma 1 (Optimal noise matrix). If Lu opt is an op- 
timal local GCP map which gives the maximum of the 
fidelity T max for the teleportation of a coherent state, 
then the resulting noise matrix is a multiple of the iden- 
tity, that is N opt = 2n op tI ■ 

Proof. First of all we observe that for a 2 x 2, symmet- 
ric and positive semidefinite matrix like N, the condition 
N = 2n opt I is equivalent to TrA = 2\/det N. There- 
fore we have to show that to pt is such that TrA opt = 
2^J det N opt . We do this by reductio ad absurdum sup- 
posing that Lu pt gives a noise matrix with TrN opt > 
2-^/det N opt . However, within the class of local GCP 
maps, there exists a subclass of local symplectic (i.e., 
unitary Gaussian) maps realized by a generic symplec- 
tic Sb on Bob mode, and the associated symplectic map 
S a = ZSbZ on Alice mode, which act as an effective sym- 
plectic transformation on N opt , N' opt = SbNoptS^ (see 
Eq. ([5|)). We can always choose Sb such that N' opt = 
y/det Nop t I, for which TrA^ pt = 2 y /detN opt < TrA opt 
while det N' opt = det N opt . However, we see from Eq. ©, 



that this local symplectic operation increases the tele- 
portation fidelity, but this is absurd because we assumed 
from the beginning that N op t is optimal. □ 

From this lemma and Eq. ^ we can therefore rewrite 
the optimal fidelity of teleportation in terms of the single 
positive parameter n opt = ^/dct N opt /2 as 



•P opt — 



1 



^opt 



(7) 



We can now derive an upper bound for T op t for a given 
entanglement of the state shared by Alice and Bob. We 
quantify such entanglement in terms of the lowest par- 
tially transposed (PT) symplectic eigenvalue, v. Such a 
parameter cannot be improved (i.e., decreased) by local 
operations and therefore provides a quantitative charac- 
terization of CV entanglement 14| . It is a local symplec- 
tic invariant and it can be expressed in terms of the four 
local symplectic invariants det A, det B, det C and det V 
as 



- 9-1/2 



v = 2 



E(V) - (S(y) 2 -4dety) 



1/2' 



1/2 



(8) 



where E(V) = det A + det B - 2 det C. 



Theorem 1 (Upper bound). For a given Gaussian 
bipartite state shared by Alice and Bob, with lowest PT 
symplectic eigenvalue v , the fidelity of the teleportation 
of a coherent state is limited from above by 



Fopt < 



1 



1 



(9) 



Proof. Let us suppose that we can achieve a larger fidelity 
T = l/(l+n op t) with < n opt < v. Alice can in principle 
have at her disposal a two-mode squeezed state, with the 
usual correlation matrix 



W = 



I cosh r — Z sinh r 
— Zsinhr / cosh r 



(10) 



(r is the squeezing parameter) and use this two-mode 
squeezed state, together with the bipartite state shared 
with Bob already optimized over all local GCP maps, to 
implement a CV entanglement swapping protocol [19| . 
In fact, by mixing at a balanced beam splitter her mode 
of the bipartite state shared with Bob and one part of 
the two-mode squeezed state, and performing homodyne 
measurements at the output, Bob mode gets entangled 
with the remaining part of the two-mode squeezed state 
in Alice hands. Since the noise added to the teleported 
state is N opt = 2n op tI, it is straightforward to see that 
the two remaining modes are then described by the fol- 
lowing CM 



W, 



swap 



I cosh r 



-Zsinhr 



-Z sinh r I[2n opt + cosh r] 



(11) 



In other words, before entanglement swapping, Alice and 
Bob shared an entangled state with CM V and entangle- 
ment characterized by v; after entanglement swapping, 
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they share a state with CM W swap . In the limit of in- 
finite squeezing the lowest PT symplectic eigenvalue of 



W, 



s<rnp tends to n op t, i.e., lim r — >00 inswap — ^opt- 
we supposed n opt < v, this means that for a sufficiently 
large squeezing parameter r, v swap < v, i.e., Alice and 
Bob have increased their entanglement. However this is 
impossible because we have employed only local opera- 
tions [20(1 . Therefore it must be n op t >v. □ 

We complete the characterization of the optimal fi- 
delity of teleportation in terms of the entanglement 
shared by the two distant parties by providing also a 
lower bound for T op t, proving in this way the result of 
Eq. ©. 

Theorem 2 (Lower bound). For a given Gaussian bi- 
partite state shared by Alice and Bob, with lowest PT 
symplectic eigenvalue v , the fidelity of the teleportation 
of a coherent state is limited from below by 



Since 



J~ op t > 



1 



1 + 3^ 



(12) 



Proof. From the definition of symplectic eigenvalue, one 
has that a 4 x 4 symplectic matrix 5* exists which diago- 
nalizes AVA (A = diag(Z, /)), i.e., the PT matrix of the 
CM V. This means SAVAS T — diag(^, v, fj,, //), where /j, 
is the largest PT symplectic eigenvalue. By writing S in 
2x2 block form 



S = 



W a W b 
W c W d 



(13) 



and rewriting the diagonalization condition for the upper 
2x2 block only, one gets the following condition 



W a Z AZWl + W a ZCW^ ' + W b C T ZWl + W b BW b T = vl. 

(14) 

The symplectic transformation S transforms the vector 
of quadratures £ into £' = (x' a ,y' a , x' b , y' b ) T = S£, and the 
PT vector into £" = {x" a ,yl,xl,y' h r ) T = One has 

[i' a ,^] = i, because commutation relation are preserved 
by S, implying 



detW a + detW b = 1. 



(15) 



The commutation relation is instead not preserved for 
the PT transformed quadratures, and introducing a real 
parameter e such that = ie, we get another con- 

dition for the two upper blocks of S, 

- detW a +detW b = e, (16) 

which together with Eq. (|15[) . gives the parametrization 

detW a = (1 - e)/2, det W b = (1 + e)/2. (17) 

Now, since Aa;" 2 = Ay" 2 = v/2, the Heisenberg un- 
certainty principle imposes that |e| < v and in partic- 
ular for every entangled state we have |e| < v < 1. 
This latter condition, together with Eq. (JT7J) , suggests 



an alternative parametrization in terms of the angle 
9 = arctan -e)/(l + e) (0 < 6 < tt/2), 



\J det W a — sin 8, \/detW b — cos 9. 



(18) 



The 2x2 matrices W a and W b and the parameter 9 
allow to construct an appropriate local map which will 
lead us to derive a lower bound for the fidelity. This local 
map is a TGCP map which, at the level of CM, acts as 



V ~>V = SVS T + G, 
with 5* and G satisfying 

G + iJ- iSJS T > 0. 



(19) 



(20) 



If the TGCP map is local, then S = S a © S b and G = 
G a © G b , with Gk+iJ- iS k JS% > (fc = a, b). 

The desired local TGCP map ue is defined in terms of 
S a , S b , G a and G b in the following way 



Sa = 

s b = 

G a 

G b = 



ZW a Z[cos9]~ L 0<6»<tt/4 . , 
ZW a Z [sm8}~ 1 n/4:<9 < tt/2, ^ & ' 



W b [cos 6] 1 0<6»<tt/4 
W b [sin^]" 1 tt/4 < 9 < tt/2, 

J [l- tan 2 9] I 0<6»<tt/4 
tt/4 < 9 < tt/2, 

< 9 < tt/4 

[1 - cot 2 9] I Tr/4<9< tt/2. [Zm> 



(21b) 
(21c) 



By applying Eqs. (fT4| and (|T9|l . one can see that this 
local TGCP map transforms the noise matrix N into a 
final matrix proportional to the identity, given by 

N = [v/cos 2 9 + 1 -tan 2 6>]J, < 9 < tt/4, (22) 
N = [tVsin 2 6> + l-cot 2 6>]/, tt/4 < 9 < tt/2. (23) 

It is however convenient to come back to the parametriza- 
tion in terms of e, which allows to express the final A in a 
unique way, for < 9 < tt/2. In fact, from Eqs. (|2"2" ]) -(|2"5 ]l . 
one gets 



A = 2- 



l + |e| 

which, inserted into Eq. ([7]), yields 



l + z/ + 2e 



(24) 



(25) 



From the condition imposed by the Heisenberg uncer- 
tainty principle < |e| < v, we see that the fidelity 
is minimum when |e| = v, so that we get the following 
lower bound 



J~ op t > 



1 + 1/ 
l + 3v' 



(26) 
□ 
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Theorems 1 and 2 provide a very useful characteriza- 
tion of the optimal fidelity which can be achieved with 
Gaussian local operations at Alice and Bob site. In fact, 
the bounds are quite tight because the region between 
the upper and the lower bound is quite small (see Fig.[T]). 
Therefore, by simply computing the lowest PT symplec- 
tic eigenvalue of the CM of the shared state and using 
the bounds, one gets a good estimate of the maximum 
fidelity that can be obtained with appropriate local oper- 
ations. In fact, the error provided by the bounds is never 
larger than 0.086 (see Fig.©. 

F 




0.2 0.4 0.6 0.8 v 1.0 



FIG. 1: (Color online) Plot of the upper and lower bounds 
(Eq. ((9} and (|12[) respectively) for the fidelity of teleportation 
of coherent states. The blue region is the allowed region in 
the (J-, v) plane. 
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FIG. 2: Plot of the distance between the upper and lower 
bounds for the teleportation fidelity versus the allowed values 
of v for an entangled state between Alice and Bob. We see 
that the error £\T under which we can estimate the maximum 
fidelity is less than 0.086. 

Corollary 1 (Upper bound achieved in the sym- 
metric case). The upper bound J- op t = 1/(1 + 1') is 
achieved iff the bipartite Gaussian state shared by Alice 
and Bob is symmetric. The optimal local transformation 
in the symmetric case is a local symplectic map. 



Proof. The "if" part of the theorem directly follows as a 
special case of the preceding proof. If the Gaussian state 
shared by Alice and Bob is symmetric, it is det W a = 
detWh, implying e = 0. Then, Eq. (j2"5| shows that in 
this case the fidelity reaches the upper bound, T op t — 
1/(1 + v). Moreover in this case 9 — ir/4 and the local 
TGCP map of Eqs. (J^J) is optimal and it is a symplectic 
one, with S a = ZW a ZV2, S b = W b V2, G a = G b = 0. 
The "only if" part instead can be easily proved by using 
the result of Theorem 3 about the CM of the optimized 
bipartite state shown in the following section. The proof 
is given in the Appendix. □ 

This latter corollary provides the generalization of the 
result of Ref. @, which obtained the same relation be- 
tween optimal fidelity and v but by considering only a 
special class of symmetric Gaussian bipartite state for 
Alice and Bob, obtained by mixing at a beam splitter 
two single-mode thermal squeezed states. 



IV. DETERMINATION OF THE OPTIMAL 
LOCAL MAP 

We have derived a lower bound for the optimal fidelity 
of teleportation of coherent states, by explicitly con- 
structing the family of local TGCP maps lob of Eq. (|21| , 
which transform Alice and Bob shared state so that 
the corresponding fidelity of teleportation is given by 
Eq. (|25|) . interpolating between the lower and upper 
bound of Theorems 1 and 2 by varying e = cos 29. The 
map tog is symplectic only for 9 — 7r/4 and in this case 
it is the local optimal map for a symmetric shared state, 
since it reaches the upper bound of Theorem 1 (see Corol- 
lary 1). When 9 ^ 7r/4, loq is a noisy (i.e., non unitary) 
map, and it is not optimal in general, because one can- 
not exclude that different LOCC strategies by Alice and 
Bob may yield a better noise matrix N and therefore a 
larger value of the teleportation fidelity. As discussed in 
the introduction, here we shall study the optimization of 
the teleportation by restricting to GCP maps, which pre- 
serves the Gaussian nature of the bipartite state initially 
shared by Alice and Bob. 

In this section we shall derive two results: i) the general 
form of the final CM of the bipartite Gaussian state after 
the optimization over all local GCP maps; ii) the opti- 
mal local TGCP map, i.e., the local TGCP map which 
maximizes the teleportation fidelity when one restricts to 
TGCP maps only, excluding in this way measurements 
of ancillary modes. 

Ref. jl3| already provided the analytical procedure for 
the determination of all the parameters of the optimal 
TGCP map. Here, by further elaborating the approach 
of Ref. [13[ , we will show that the optimal TGCP map 
can always be written in a simple form, as a local sym- 
plectic operation eventually followed by a single mode 
attenuation [26[, either at Alice or Bob site. 
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A. Standard form of the correlation matrix of the 
optimized bipartite state 

In this subsection we show that, even if we do not know 
the specific form of the optimal GCP map, one can always 
characterize it indirectly by determining the general form 
of its outcome, i.e., the general form of the CM of the 
final Gaussian state shared by Alice and Bob after the 
maximization. We begin with the following lemma. 

Lemma 2 (Standard form III). The correlation ma- 
trix V of every bipartite Gaussian state can be trans- 
formed by local symplectic operations into the following 
normal form 



Vl 




(27) 



where all the coefficients are positive and satisfy the fol- 
lowing constraints 



m - n 2 
That is: 



mi — m 2 = di 



n + X 



A. 



-d- A 



A e 



Vi = 



-d — A m + A 



(28) 



(29) 



Proof. It is well known that it is possible to transform 
every V of an entangled state in the usual normal form 
(standard form I) [23, H| 



V N = 



\ 



-ci 



C2 




(30) 



where all the coefficients are positive. Now we per- 
form a local symplectic operation composed by two lo- 
cal squeezing operations, S a = diag(y / r^", 1/y/rZ) and 
Sb = di&g(^/r b , l/y/r~b). We impose the first two con- 
ditions n\ — n 2 — mi — mi = A, 



"K-^ 1 ) = Kn-r b l ) = A, 

which solved for positive r a and r b give 



(31) 



f a (A) = A/(2o) + y/1 + (A/2a) 2 



r b (X) = X/(2b) + + (A/26) 2 . 
Now we impose the last constraint d\ — d 2 = 

ci^r a {\)r b {\) - c 2 /y/r a (X)r b (X) = A. (34) 



(32) 
(33) 

A, that is 



Our lemma is proved if there is at least one solution A of 
Eq. (|3"4"|) . Since A = is the trivial solution when ci = c 2 
and Vjv = Vi, we can exclude this particular case and 



divide Eq. (|34|) by A. Therefore we have to show that the 
equation 

1 



/(A) = j[ci^/r a (X)r b (X) - c 2 / ^ r a {X)r b {X)} = 1 (35) 

admits at least one real solution. If |A| 3> a and |A| ^> 6, 
we can power expand the square roots in Eqs. (I32p ~ (l33p 
so that we easily find the following limits for /(A): 



lim /(A) 

A — *OC 



ci/Vab < 1, 
lim /(A) = c 2 /Vab < 1, 

±sign(ci - c 2 )oo. 



lim /(A) 



(36) 
(37) 
(38) 



The inequalities in Eqs. (|36 p - ([37|) follow from the Cauchy- 
Schwartz inequality (xi)(x 2 ) > (xix 2 ) 2 applied to the 
quadrature operators of the two modes. Given the three 
limits (|36|) . ([37f and (|38|) . since /(A) is continuous ev- 
erywhere except at the origin, at least one solution of 
Eq. (|55|) exists. Moreover this solution has the same sign 
of ci — c 2 . □ 

We have defined the standard form of lemma 2 as stan- 
dard form III because it is very similar to the standard 
form II defined in Ref. [27| for the determination of a nec- 
essary and sufficient entanglement criterion for bipartite 
Gaussian states. In particular the two standard forms co- 
incide in the special case of a symmetric bipartite state 
(ni — mi and n 2 — m 2 or equivalently n = m). 

Theorem 3 (Form of the CM of the optimized 
bipartite state). The optimal GCP map Lo op t maximiz- 
ing the teleportation fidelity is such that the CM of the 
transformed bipartite state is in the standard form III V± 
defined by Eqs. @7j)-f29j). 

Proof. By means of local symplectic operations, we can 
always put the CM of the bipartite state of Alice and Bob 
in the form of Eq. (|2"7| , but without the constraints of 
Eq. I|28p. We first restrict to local symplectic operations 
and show that the optimal local symplectic operation al- 
ways transforms to a state with a CM satisfying the con- 
straints of Eq. (|28p. Since the CM is tridiagonal, then 
any possible optimal map must be a squeezing transfor- 
mations of the two modes given by S a = diag(r a ,r~ 1 ) 
and Sb = diag(r&, r^ 1 ) [l3j]. Let us define 



a(r a ,r b ) 
P(r a ,r b ) 



r 2 a ni 



2dir a r b + r b mi, (39) 
r- 2 n 2 - 2d 2 (r a r b )- 1 + r b 2 m 2 , (40) 



so that the noise matrix of Eq. is equal to N — 
diag(a, (3). The optimal map must minimize det N = a/3, 
and therefore we impose that 



Va(r a ,r b )f3(r a ,r b ) = 0, 



(41) 



where V = {d ral d rb ). Due to Lemma 1, we must 
also have that a(r a ,r b ) = f3(r a ,r b ) ^ 0, and therefore 
Eq. (|4T"T) reduces to 



V[a{r a ,r b ) + f3(r a ,r b )] = 0. 



(42) 



7 



The CM of the transformed state is the optimized one 
iff the optimal local symplectic operation is the identity 
map, that is, if a(l, 1) = (3(1, 1) and 



V[a(r a ,r b ) + (3(r a ,r b )} = 0. 

r a =r h = l 



(43) 



It is easy to check that these conditions are satisfied iff 



Til — 77.2 = 777-1 — 7712 



di 



(44) 



which are exactly the constraints of Eq. The the- 

orem is proved if we show that the normal form V\ of 
Eq. is actually kept also if one maximizes over the 
broader class of GCP maps. In fact, if by reductio ad 
absurdum, we assume that an optimal (non-symplectic) 
GCP map exists leading to a CM not satisfying the con- 
straints of Eq. (f2"8")) . we could always apply a further 
symplectic map which, by repeating the maximization 
above, would transform to a state with a CM satisfying 
the constraints ([28| and yielding a larger teleportation 
fidelity. But this is impossible, because it contradicts the 
initial assumption of starting from the optimal bipartite 
state. □ 

In other words, the CM of the state shared by Alice and 
Bob after the maximization of the teleportation fidelity 
is the one with the standard form III of Eq. because 
it is the unique CM for which the optimal map is the 
identity operation on both Alice and Bob site. 



B. Optimal trace-preserving Gaussian CP map 

In the former subsection we have determined the form 
of the CM of the optimized state of Alice and Bob, with- 
out determining however which is the local GCP map 
which maximizes the teleportation fidelity. Here we find 
this optimal map, restricting however to the smaller class 
of trace-preserving GCP maps. The case of Gaussian 
maps including Gaussian measurements on ancillas will 
be afforded elsewhere. 

Ref. [221 ] has introduced the notion of minimal noise 
TCGP maps, as the extremal solution of the condition of 
Eq. (|20[) . These maps are the ones that, for a given ma- 
trix S, possess the "smallest" positive matrix G realizing 
a CP map. It is easy to check that a minimal noise TGCP 
map satisfies the relation det G = (1 — det S) 2 . An exam- 
ple of minimal noise TGCP map is an attenuation [26| . 
i.e., the transmission of a single boson mode through a 
beam splitter with transmissivity r (0 < r < 1), such 
that 



T 2 a v , 



(45) 



where a = (x + ip) /y/2 is the annihilation operator of 
the mode, and ay that of the vacuum mode entering the 
unused port of the beam splitter. 

It is evident that the TGCP map maximizing the tele- 
portation fidelity has to be a minimal noise TGCP map 
131. We prove now a useful decomposition theorem. 



Theorem 4 (Decomposition of TGCP maps). A 

minimal noise TGCP map to on a single mode system 
with det G < 1 can always be decomposed into a symplec- 
tic transformation o\, followed by an attenuation t and 
by a second symplectic transformation o~2, that is, 



CO 



0~2 O T O 0~\. 



(46) 



Therefore a local minimal noise TGCP map on a bipartite 
CV system can always be decomposed into a local sym- 
plectic map, followed by the tensor product of two local 
attenuations and by a second local symplectic map. 

Proof. We consider a generic minimal noise TGCP map 
such that V V = SVS T + G for a generic CM V, with 
det G = (1 — det S) 2 . G is a positive symmetric matrix, 
and therefore a symplectic matrix T2 exists such that 
G = T 2 T 2 T (l - s), where s = det 5". We then define the 
symplectic matrix T\ = -^T^S, and we also consider 

an attenuation map with transmissivity y^s. If we now 
first apply the symplectic map defined by T\, then the 
attenuation map and finally the second symplectic map 
defined by T2, by using the relations Ti\fsT\ = S and 
G = T^T/^l — s), one can check that the composition of 
the three maps reproduces the given TGCP map. □ 

Corollary 2. A minimal noise TGCP map to on a single 
mode system with G proportional to the identity matrix, 
i.e., G = (1 — s)I (0 < s < I) can always be decomposed 
into a symplectic transformation o\, followed by an at- 
tenuation T, 



T O 0~\. 



(47) 



Therefore a local minimal noise TGCP map on a bipartite 
CV system with Gi = (1 — Si)I (0 < Sj < 1, i = a, b) 
can always be decomposed into a local symplectic map, 
followed by the tensor product of two local attenuations. 

Proof. It is sufficient to repeat the former proof and con- 
sider that, since G = (1 — s)I, T2 = I and therefore the 
second symplectic map is the identity operation. □ 

This latter case is of interest because the optimal 
TGCP map must have in fact the property Gi — (1 — Si)I, 
i = a,b. To show this we first s imp lify the scenario 
by exploiting the results of Ref. [13j . which provides 
the general analytical procedure to derive the optimal 
local TGCP map. In fact, Ref. [HI shows that when 
the optimal local TGCP map is a minimal noise, non- 
symplectic one, it can be performed on one site only, 
i.e., either on Alice or on Bob alone. Suppose that 
the non-symplectic map is performed on Bob site; it is 
straightforward to see that, under a generic local TGCP 
map V -> V = SVS T + G, with S = I © S b and 
G = I G b , the noise matrix N transforms according 
to N -> N' = T + G b where 

V = ZAZ + ZCSj + S b C T Z + SbBSj . (48) 
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Ref. [H[ shows that the optimal map is such that T oc 
Gf,, and since Lemma 1 shows that the optimal N is 
proportional to the identity, this implies that Gb must 
be proportional to the identity. Therefore Corollary 2 
leads us to conclude that the optimal local TGCP map 
is either a local symplectic map, or a local symplectic 
map followed by an attenuation by a beam splitter, placed 
either on Alice or on Bob mode. 

Theorem 4 and Corollary 2 therefore provide a very 
simple and clear description of the TGCP map which 
maximizes the teleportation fidelity, which is not evident 
in the treatment of Ref. [l3| . We can further characterize 
the optimal local TGCP map by determining: i) the form 
of the first local symplectic map; ii) the conditions under 
which the optimal local operation is noisy, i.e., when one 
has also to add a beam splitter with appropriate trans- 
missivity on Alice or Bob mode in order to maximize the 
teleportation fidelity. In order to do that we first need a 
further lemma, similar to lemma 2. 

Lemma 3. For any given positive real parameter rj, the 
correlation matrix V of every bipartite Gaussian state 
can be transformed by local symplectic operations into the 
following normal form 



Yn 




(49) 



where all the coefficients are positive and satisfy the fol- 
lowing constraint 



n 1 -n 2 = r](di - d 2 ) = rj 2 (m 1 - m 2 ) = A, 



A e 



(50) 

That is, we have a family of normal forms depending on 
the parameter rj, 



V n = 



I n + A — d — A/77 

n d 
—d — A/77 m + ^/v 2 

\ d 771 



(51) 



Proof. With the same procedure used in the proof of 
Lemma 2, we arrive to an equation similar to (|35|) . while 
the corresponding three limits are exactly the same of 
(|3"T|) and ([55]) . since the factor 77 cancels out. As a 
consequence the continuity argument is valid also in this 
case, and therefore, for every fixed parameter 77, one can 
find a transformation which puts the CM in the normal 
form V v . □ 

We notice two facts that will be useful for the next 
theorem: i) the optimal CM standard form of Theorem 
3, Vi, belongs to the class of normal forms V n , since it is 
obtained for 77 = 1; ii) when < 77 < 1, the state with 
CM V v is transformed into the Gaussian state with CM 
equal to V\ when a beam splitter with transmissivity 77 is 
put on Bob mode. We then arrive at the theorem about 
the optimal TGCP map. 



Theorem 5. The optimal local TGCP map maximizing 
the fidelity of teleportation of coherent states can always 
be decomposed into a local symplectic map, eventually fol- 
lowed by an attenuation either on Alice or on Bob mode. 
The first local symplectic map is the one transforming 
the CM of the Gaussian state shared by Alice and Bob 
into one particular normal form of the family V n defined 
in Eqs. ^49\ )- [oU\) , with < 77 < 1. One has to add the 
attenuation on one of the two modes for realizing the op- 
timal TGCP map if there is a value oj "rj, let us say 77 = t. 
such that the coefficients of V T satisfy the relations 



d 



1' 



< 1. 



(52) 



// instead the condition of Eq. I52j) is never satisfied in 
the interval r\ G [0, 1), the optimal TGCP map is formed 
only by the local symplectic map transforming the CM 
into the normal form Vi (i.e., 77 = 1 and no attenuation 
is required). 

Proof. Using Lemma 3, we can always apply a local sym- 
plectic map which transform the CM into one of the form 
of the family of Eq. (|5ip . with 77 = r. We then apply an 
attenuation on Bob mode with transmissivity r and then 
try to find the maximum fidelity as in the proof of The- 
orem 3. We have that the two diagonal elements of the 
noise matrix N now read 

a(r) = /3(t) =n-2Td + T 2 m+l-T 2 , (53) 

The optimal map must minimize dct N = a 2 , and there- 
fore we impose 



da(r) 
dr 



= 0, 



(54) 



which is satisfied iff t = dj (m— 1). If < t < 1, this map 
composed by the local symplectic map and the attenua- 
tion is the optimal map. If instead for any r G [0, 1], the 
condition of Eq. (p)2")) is not satisfied, there is no critical 
point in this interval and therefore the optimal map is 
just the symplectic transformation to the normal form 
V\. One has also to check the behavior at the lower 
boundary value r = 0, but this is trivial because this 
means that Bob uses the vacuum to implement the tele- 
portation which is never the optimal solution if we have 
an entangled channel. In fact, if Bob uses the vacuum, 
the channel looses its quantum nature and the maximum 
of the fidelity is the classical one T — 1/2, which is be- 
low the lower bound for any entangled state given by 
Eq. d). □ 



Theorem 5 therefore characterizes in detail the optimal 
TGCP map, giving in particular the conditions under 
which this map is noisy, i.e., non-symplectic and therefore 
when teleportation is improved by increasing the noise 
and decreasing the entanglement of the shared state. 

Using this latter theorem we can also determine how, 
from an operational point of view, one can compute 
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the value of the tclcportation fidelity maximized over all 
TGCP maps, starting from the symplectic invariants of 
the bipartite Gaussian state initially shared by Alice and 
Bob. From the CM of this latter state one can: 

1. compute the four symplectic invariants a = 
Vdet A, b = Vdet B, c = y/\ det C\ and v = det V. 

2. Knowing the first three invariants of the channel, 
the elements n, m, and d of the normal form V n can 
be expressed as functions of only the two unknown 
parameters A and r\ as 



n(A) = -A/2 + v/« 2 + (A/2) 2 , 
m(A,r?) = -\/2T? + y/iP 



(X/2rj 2 



-X/2ri + \/c 2 + (A/2??) 2 



(55) 
(56) 
(57) 



3. The two parameters A and rj can be found solving 
the following system in the region < rj < 1, 



detV n (X,r]) = v 
V [m(A, 77) - 1] = d(X,t]) 



(58) 



and solving also the first equation in the boundary 

T) = 1, 



det Vi(A) = v, 



(59) 



(the two conditions of ([55)) come from the invari- 
ance property of the determinant of the channel 
and form the maximization condition of Eq. I|52p ). 

4. We call (A^, rji) with i — 1,2 ... k, the union of the 
solutions of (|5"5|) and (JHSJ) (we have at least one 
solution because (fB"9"]) admits at least a solution). 
We then compute the candidate fidelities 



Ft = 2 



2 + V / « 2 + A f/4+^f + A2/4 



2J<?r,l + Xl/4+l-r,i 



(60) 



so that the maximum fidelity will be F op t 
max{Fj}. 



V. CONCLUSIONS 

We have studied the Braunstein-Kimble protocol Q 
for the CV teleportation of coherent states and how the 
corresponding fidelity can be maximized over local oper- 
ations at Alice and Bob site. We have assumed that Alice 
and Bob share a Gaussian bipartite state and restricted 
to Gaussian LOCC strategies, which preserves such a 
Gaussian property. We have shown that, for a given 
shared entanglement, the maximum fidelity of telepor- 
tation is bounded below and above by simple expressions 
depending upon the lowest PT symplectic eigenvalue v 



only (see Eq. iJTJ). We have seen that these bounds are 
quite tight and that the upper bound of the fidelity is 
reached if and only if Alice and Bob share a symmet- 
ric entangled state. We have also determined the general 
form of the CM of Alice and Bob state after the optimiza- 
tion procedure. Then we have restricted to local TGCP 
maps and shown that the optimal TGCP map is com- 
posed by a local symplectic map, eventually followed by 
an attenuation either on Alice or on Bob mode. Finally 
we have shown how the corresponding value of the max- 
imum fidelity F op t can be derived from the knowledge 
of the symplectic invariants of the initial CV entangled 
state shared by Alice and Bob. 

If one considers generic GCP maps (i.e., including also 
Gaussian measurements on ancillas), one expects to fur- 
ther improve the teleportation; in this case one should 
adopt the description of GCP maps given in [24| in order 
to characterize the optimal local Gaussian map, but this 
will be the subject of future work. 

Another open question is to see if teleportation fidelity 
can be increased by leaving the Gaussian setting stud- 
ied here and consider more general non-Gaussian local 
operations at Alice and Bob site. In this respect, the 
preliminary results of Refs. [25j seem promising. 
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VII. APPENDIX 

We now prove the "only if" part of Corollary 1, i.e., 
that if J- opt = 1/(1 + 1/) (the upper bound of the opti- 
mal fidelity) , then the bipartite Gaussian state shared by 
Alice and Bob is symmetric. 

Proof. Theorem 3 shows that the final CM after any op- 
timization map must be V\ of Eq. (|29p . Using Lemma 1 
and the explicit form of V\ , the hypothesis is equivalent 
to n + m — 2d — 2v 1 so we can make the substitution 
d = (n + m)/2 — v, which is just a different parametriza- 
tion: V\(n,m,d,X) — > V\ (n, m, f, A) . Now, the condition 
that v is equal to the PT minimum symplectic eigenvalue 
gives us a constraint on the parameters of the matrix V\ 



u{Vi(n, m, v, A)} = v. 



(61) 



If the state is symmetric, which means n = m, then the 
condition of Eq. (f6T|) is identically satisfied. If the state 
is non symmetric, Eq. (|61[) is a non-trivial equation that 
solved for A gives A = (m — n) 2 /8f — n — m. However the 
corresponding matrix V\ = (n, m, v, X) is not the CM of 
a physical state; in fact, the characteristic polynomial of 
V\ can be written as P{x) = {cQ + cix + x 2 )(gQ+ gix + x 2 ) 
where cq = —v(n + m + v), but this means that V\ has 
at least one negative eigenvalue and therefore it is not 
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positive definite. Therefore T op t — 1/(1 + v) is realized 
only if Alice and Bob state is a Gaussian symmetric state. 
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